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Some key points out of my recent study of the characteristic features of the small (x<g;l), intermediate (:r~l) 
and large (x S> 1) Leutwyler-Smilga regimes for QCD-type theories in a finite volume (x = VSm) are presented, 
and a few immediate consequences for full QCD simulations are discussed. 



be an option; hence one is left with the question 
whether the given degree of (non-)ergodicity is 
tolerable for practical purposes or whether an en- 
semble average on the basis of the data available 
is afflicted with a sizable systematic error. It is 
clear that the answer may depend on the observ- 
able one has in mind. Hence from a physical point 
of view the question may be rephrased as follows: 
Does a given observable on a perfect (full) QCD 
ensemble - generated in a box of given volume V 
and with Nj sea-quarks of mass m - depend on 
the topological charge v of the background ? 

2. Leutwyler-Smilga analysis in a nutshell 

For pionic observables, much can be said on 
purely analytical grounds. Leutwyler and Smilga 
(LS) have shown [Bj that a clear answer is possible 
for two regimes of quark masses and box volumes, 
both of which involve the LS-parameter[] 




Figure 1. Time history of zy na i in a HMC simu- 
lation with a quadruple of staggered quarks at 
(3 = 5.35, m — 0.01 (figure taken from no. 3 in jij). 

1. Motivation 

Numerical simulations of full QCD may be 
hampered by the fact that certain observables 
show rather long autocorrelation times. A well 
known example is the global topological charge 
v which was found to develop problems for 
sufficiently small quark masses: For staggered 
fermions an insufficient tunneling rate has been 
observed both in the symmetric and in the bro- 
ken phase [0. For Wilson type sea-quarks mo- 
bility between the topological sectors was found 
to be satisfying for M^/Mp > 0.56, but the same 
problem emerges if k is tuned sufficiently close 
to K cr ;t 0. Fig. 1 displays the time-history of a 
simulation representing the "worst case scenario" ; 
the sample is obviously non-ergodic, since the 
generated ^-distribution is far from symmetric. 

Simply prolonging the simulation time may not 
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which indicates whether the box is "small" or 
"large" in the sense that it decides on whether 
the systems prefers to show symmetry restora- 
tion or SSB phenomena, respectively: For x <C 1 
chiral symmetry is effectively restored and quarks 
and gluons represent appropriate degrees of free- 
dom. For x^>l the SU(Nj)a symmetry is quasi- 
broken, meaning that - though the box-volume is 
formally finite - pions represent appropriate de- 
grees of freedom (Nf > 2). Note that m refers 
to the mass of the sea-quarks and that x is (in 
principle) independent of the parameter M n L. 

1 S = — lim m ^o liniy_ hoo (ijiip) (this order) is the conden- 
sate in the chiral limit; both E and m are scheme- and 
scale-dependent, but the combination is an RG-invariant. 
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Figure 2. Distribution of v gco in the small (8x4 lattice), intermediate (18x6) and large (40x10) Leutwyler- 
Smilga regimes, respectively [QED(2) data for = 3.4, m = 0.09, Nf = 2]. 



The net result of the LS-analysis [|3| is that for 
x <C 1 the partition function is completely domi- 
nated by the contribution from the topologically 
trivial sector for small enough m, since 
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In the opposite regime i>1 and with the addi- 
tional condition that the quark masses are so light 
that the sea-pion overlaps the box (M^L^ 1), the 
path-integral in the effective description (requires 
A c gL^l) is dominated by the constant mode 
and Leutwyler and Smilga end up finding that the 
partition function is approximately independent 
of v or, more precisely, that the distribution is H 
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3. Checking the net charge distribution 

Since the LS-issue is peculiar to the full theory, 
a preliminary study in a QCD-type theory seems 
permissible. The massive multiflavour Schwinger 
model (QED(2) with Nf > 2) is supposed to repro- 
duce all qualitative features ||. I use the Wilson 
gauge action i5 gaU g e =/3$^(l — cosf? n ) and a pair 
of (dynamical) staggered fermions. 

I compare the three regimes x-^il, i~1,s>1 
to each other using three dedicated simulations: 
Working with fixed staggered mass m— 0.09 and 
(3 = 1/g 2 — 3.4, the three regimes are represented 
by the three volumes V = 8 x 4, 18 x 6, 40 x 10. 
From this setup the LS-parameter takes the val- 
ues x~0.33, 1.12,4.16, respectively, and the pion 
(pseudo-scalar iso-triplet) has a (common) mass 
M„ = 0.329 and therefore a correlation length 
^tt =3.04 as to fit into the box. 



A configuration is assigned an index only if 
the geometric (f gC o = 4z Yl 1°S an( i the field- 
theoretic definition (vith = K v a aii v n&i = Yl, sm @o > 
k ~ 1/(1 — (S g augc ) / f3V), after rounding to the 
nearest integer, agree — for details see |q|. 

As one can see from Fig. 2, the charge distri- 
butions found in the three runs seem to follow a 
general pattern consistent with the LS-prediction: 
For x <§; 1 the distribution (and therefore the par- 
tition function) is dominated by the contribution 
from the topologically trivial sector, whereas for 
i>1 the distribution gets broad and seems com- 
patible with the gaussian form (||) with variance 



c/Nf, i. 



-{v 2 )=x/2. 



4. Role of the functional determinant 

Since the LS-issue is peculiar to the full (un- 
quenched) theory, an attempt to understand by 
which mechanism the three regimes differ from 
each other may lead one to investigate how the 
functional determinant or specifically its contri- 
bution to the total action per continuum-flavour 

Sformion = ~ log(det($> + 77l)) + COnSt (4) 

relates to the contribution from the gauge field 
and, in addition, how this might depend on the 
topological charge of the background. The idea 
is thus to study such a relationship sectorally, i.e. 
after the complete sample has been separated into 
subsamples with definite absolute charge \v\. 

From the scatter plots in Fig. 3 one notices a 
(weak) positive correlation between Sformion and 



This means that the functional determi- 



nant acts - roughly - like an effective renormal- 
ization of (3 with a factor bigger than 1. A key ob- 
servation is that the correlation improves, if one 
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Figure 3. Sf e rmioii = — log(det(^+m)) (i.e. per continuum flavour) versus S gauge together with the best 
linear fits for the reasonably populated topological sectors in the small (8x4 lattice), intermediate (18x6) 
and large (40x10) Leutwyler-Smilga regimes, respectively [QED(2) data for /3 = 3.4, m= 0.09, Nf=2]. 



separates the sample into subsamples with fixed 
\v\ — for each one of which the best linear fit for 
Sfcrmion as a function of "S'gaugc is included. Quite 
generally, the linear function for low \v\ lies below 
the one(s) for higher \v\, and has a larger slope 
than the latter. This means that the functional 
determinant brings - in general - an overall sup- 
pression of higher topological sectors w.r.t. lower 
ones and a sectorally different renormalization of 
(3. There are interesting effects as x varies: In the 
small LS-regime both the sectoral dependence of 
the renormalization factor and the offset between 
neighboring sectors are huge, i.e. for x — > the 
functional determinant acts as a constraint to the 
topologically trivial sector. For x ~ 1 both offset 
and sectoral dependence of the slope are smaller. 
In the large LS-regime there is only a minor over- 
all suppression of higher topological sectors w.r.t 
lower ones, and the effective renormalization of (3 
(and hence, in 4 dimensions, of the physical lat- 
tice spacing) seems to be uniform for all sectors. 

5. Sectoral heavy-quark potential 

It is instructive to play the same game with a 
standard observable like the heavy-quark poten- 
tial. As one can see from the leftmost column 
in Fig. 4, the physical potential (on the unsepa- 
rated sample) shows, by its curvature, clear evi- 
dence of the dynamical fermions - the quenched 
potential would be strictly linear. The remain- 
ing columns show for each LS-regime the sectoral 
heavy-quark potentials on the subsets with low- 
est and highest (reasonably populated) \v\. The 
sectoral dependence is found to be striking in the 
regime i < 1, where the two subpotentials are 



clearly inconsistent with each other. For inter- 
mediate x, there is still a clear disagreement be- 
tween the two extreme- 1^| subpotentials shown. 
In the regime i>1, finally, data show good con- 
sistency between neighboring topological sectors 
(e.g. v = and v = ±l - see [||), but there is no 
clear evidence of an overall consistency. In fact, 
the sectoral potentials for v = and v = ±3 seem 
to indicate that even in the large LS-regime a 
remnant sectoral dependence might persists upon 
generating extreme changes in the absolute value 
of the topological charge. Of course it is pos- 
sible that in the large- a; simulation at hand the 
value of x (about 4.16) is not yet large enough for 
the large- a; behaviour to be fully pronounced, but 
the rule of thumb seems to be that in the large- a; 
regime standard physical observables may be sen- 
sitive to distortions of the ^-histogram which are 
large compare to its natural width cr~ (x/Nf) 1 / 2 . 

6. Lessons for full QCD simulations 

Each one of the LS-regimes has its characteris- 
tic features: For x<^l correct sampling w.r.t. v is 
easy to achieve; a simple constraint to the topo- 
logically trivial sector proves useful. For x ~ 1 the 
situation is much more delicate, since this is the 
regime where good sampling w.r.t. the topological 
charge is both crucial and nontrivial to achieve. 
For x^>l standard observables (i.e. unrelated to 
the £/(l)A-issue) prove fairly insensitive on the v- 
distribution — unless the latter suffers from ex- 
treme distortions, i.e. on a scale which is large 
compared to its natural width cr~ (x/Nf) 1 / 2 . 

In spite of the good news that state-of-the-art 
simulations of full QCD operate in the large-x 
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Figure 4. Physical heavy-quark potential (leftmost column) and selected sectoral versions in the small 
(V = 8x4), intermediate (18x6) and large (40x10) Leutwyler-Smilga regimes, respectively [QED(2) data]. 



regime j(| , there is - in principle - a problem with 
chiral extrapolations. It derives from the fact that 
the key role of the LS-parameter is to determine 
which one of the limits to — > and V — ► oo is in the 
"inner" position and thus winning: x 3> 1 corre- 
sponds to lim m ^o hmy^co , whereas a; <C 1 imple- 
ments a situation like limy^oo lim m ^o- There- 
fore one likes to extrapolate from data-points won 
from a cascade of simulations for which x in- 
creases (or stays at least constant and ^> 1) as 
the dynamical quark mass decreases. Looking at 
the definition (Q), one realizes that x decreases if 
a series of simulations on a given grid is used, i.e. 
one is marching in the wrong direction: If f3 gets 
reduced for smaller to so as to keep the lattice 
spacing constant (cf. Fig. 3 and the subsequent 
discussion), x decreases only in proportion to to; 
if one works at fixed j3, the situation is of course 
worse. To correct, one would have to use larger 
grids for smaller quark masses. While in current 
state-of-the-art simulations the box to start with 
seems sufficiently large (or the quark mass to end 



up with heavy enough) as to prevent this effect 
from getting numerically important || , it is clear 
that this is a point one may wish to control as the 
quest for lighter quark masses goes on. 
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